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Abstract 

The transition amplitude for ttK scattering is evaluated within the 
SU(3) Nambu-Jona-Lasinio model. Ordering terms according to the 

Ph. 

<-j ! expansion in 1/N C leads to a box-like diagram, t channel diagrams that 

admit scalar isoscalar (a, a') exchanges, and a u channel exchange of a 

X. 

scalar isodoublet ax that has quantum numbers corresponding to the 
Kq(1430). Both the Pauli-Villars and 0(3) regularization procedures 
are used to evaluate the T = values of the I = scattering lengths 
Oq 7 ' 2 and clq^ 2 . The finite temperature dependence is studied. We find 
that the variation in the t channel in the calculation of a^ 2 leads to a 

3 /2 

change in a Q of a factor of about two over the temperature range of 
T = 150 MeV. 
PACS: 11.30.Rd; 13.75.Lb; 14.40Aq; ll.10.Wx. 
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I. INTRODUCTION 



The pseudoscalar octet (jr, K, 77) forms the set of lowest mass hadronic states and 
are thus particles that play a central role in heavy ion collisions, since they occur 
copiously In particular the pion, with lowest mass and being an isospin triplet 
has a central role, while the kaonic sector, occurring as isospin doublets, gives us 
information on strangeness production. Since this octet consists of the Goldstone 
bosons associated with spontaneous chiral symmetry breaking, their role in nature 
is determined almost exclusively by this feature. Conversely, our understanding of 
chiral symmetry breaking is augmented by studying the properties of such mesons, 
and this can be extended to finite temperatures and baryon density. 

The purpose of this paper is to examine elastic irK — > 7rK scattering in particular 
at T = and at finite temperatures. Expressions for the transition amplitudes 
and hence the scattering lengths at T = are derived and then the temperature 
dependence of the scattering lengths is investigated. This is of interest (a) in its 
own right, although at present it is unclear how one would access finite temperature 
scattering lengths experimentally, and (b) as it provides a first step in the calculation 
of the elastic tcK — > ttK cross-sections as a function of the temperature that are 
required as input for the collision dynamics for a simulation of heavy-ion collisions 
using a chiral Lagrangian. 

At T = 0, ttK scattering has been studied within the framework of chiral per- 
turbation theory (CHPT) ]IJ-0|. It is not a priori clear to what extent a chiral 
expansion is appropriate for properties involving the strange quark, in view of its 
mass. Yet the results of even the early calculations of the scattering lengths^] ||T| 
fall within the expected measured range |)]-[7j , although one should note that actual 



There is an early Russian publication on this subject that we also draw the reader's 
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experimental values are rather poorly known. An elaboration of this approach that 
includes intermediate resonances || , in particular the exchange of the vector K* in 
the w-channel as well as of a p-meson in the t-channel, alters the actual numerical 
values by 2% at most. 

In this paper, we wish to examine the ttK scattering lengths at finite temperature. 
The Nambu-Jona-Lasinio (NJL) model gives a good description of these deeply 
bound mesons at T = 0. In obtaining such a description, the mean field or Hartree 
approximation to the self-energy, together with the random phase approximation for 
the mesonic fields has been used. It has been recognized that these approximations 
taken together constitute an expansion in the inverse number of colors 1/N C , within 
the model 0,0] and are essential in order to preserve the chiral symmetry of the 
underlying Lagrangian. For T / 0, the dependence of some of the pseudoscalar 
and scalar mesons as a function of temperature has qualitatively been validated 
by lattice gauge calculations [I2]jl3|] . The results of Ref. JL3J display a remarkable 



qualitative agreement with the finite temperature behavior of the NJL model (pole) 
masses. This gives one some confidence in the application of this model in this form 
to finite temperature, although it is speculative. 

In addition, it can be shown that the current algebra results or alternatively 
the results of chiral perturbation theory to lowest order can be obtained from this 



model on making a suitable chiral expansion [II. K|. It is then a simple issue to 
extend calculations to finite temperature, and examine the consequences thereof. 
There is no necessity to enforce the chiral limit, since this is a model study. We thus 
examine the three flavor NJL model. Contributions to the scattering amplitude are 
organized according to the expansion in 1/N C . To lowest order within this scheme, 
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there are several sorts of graphs that occur: a box graph is present, analogous to 
the four-point tree-like structure of mesonic theories. In addition, in the t channel, 
exchange of neutral scalars, in this case the <Jq and a$ — > o and a' states occur. 
This has no counterpart in the CHPT version |. Finally, the u channel contains 
the exchange of a ax, corresponding to the Kq meson - note that this also differs 
from the CHPT calculations that explicitly include vector resonance exchange of the 
K* [0J. Note that in our case, no six point tadpole-like diagrams are present. In 
addition, unitarising graphs that include irK scattering in the intermediate state are 
not included to this lowest order in 1/N C calculation of the scattering lengths in the 
NJL model, but which would be essential for a complete description of the scattering 
amplitude of cross section calculated as a function of center of mass energy. 

Before discussing our results, we comment on the procedure that is to be taken, 
in particular with regard to the application of the mean field plus random phase 
approximation at T = and at T ^ 0. As discussed earlier, this coupled symmetry 
preserving approximation works excellently at T = 0. This is true for calculating 
static mesonic properties, such as masses, mean charge radii and scattering lengths, 
particularly in the SU(2) sector, where these have been studied in detail and are 
well-known. Corrections to the mean field values introduce chiral logarithms, see for 
example [IB] and further powers in the pion mass. Calculations of the full scattering 



cross-sections, however, do require these corrections, even at T = 0, in order to 
preserve unitarity. Thus only scattering lengths can be reliably calculated. This is 



2 Note that no p-meson can be exchanged in this minimal version of the NJL model, but 
would also be of the same order as the a exchanges in the 1/N C expansion within an 



extended NJL model [15]. For this reason, we do not calculate the I = 1 scattering lengths 

1/2 3/2 
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well-established for the two flavor case [fL7[ and is the reason for our restriction to 
investigating this quantity only. 

In extending our calculation to finite temperature, we are being speculative, 
encouraged however, by the success of the qualitative agreement of the lattice pre- 
dictions of the temperature variation of particle masses with those predicted by 
the model. Furthermore, it is not expected that terms beyond this leading 1/N C 
calculation play a dominant role, since the chiral logarithms and powerlike mass 
terms that might then occur depend then on the temperature dependent pion mass, 
which in turn stays constant within the model for a wide range of temperatures, 
even up to 200 MeV. The divergences that have been found in this work in the 
scattering lengths are induced by the presence of the inverse of the decay constants 
which themselves become small at high temperatures. This feature should remain 
unchanged by higher order corrections. 

Our model results depend on the regularization procedure used, and can therefore 
at most be qualitative. While the Pauli-Villars method is generally preferred at 
T = since it preserves all symmetries, and enables one to check that the crossing 
symmetries are correctly implemented, at T / 0, Lorentz invariance is broken by 
the heat bath so that a three dimensional cutoff can be used. We thus only expect 
to obtain a qualitative description. We evaluate the scattering lengths at T = 
using both prescriptions, and considering parameters sets of other authors |19| , p0| ,p[] 
that were used for different purposes to fix other variables. In this sense, we do not 
fit model parameters so as to explicitly give us good values of the scattering lengths. 
Our results at T = are obtained using both regularization schemes. At finite 
temperatures, we find that the t channel varies somewhat with temperature due 
to the exchange of the scalar mesons, which themselves have a strong temperature 
dependence, and we find that the variation of a^ 2 in particular, is visible. At 
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T = 150 MeV, the T = value is increased by a factor of 2.36. This is in contrast 



to the temperature dependence found for the irn scattering lengths |2l| which show 

1/2 

almost no change over T = to 150 MeV. The value of a changes only slightly, 
increasing in this case only by a factor of 1.18. 

This paper is organized as follows. In Section ||, we give a brief introduction to 
the NJL Lagrangian, sufficient to introduce our notation. In Section PH we classify 
the scattering diagrams required for elastic ttK scattering within the framework of 
the model and examine in particular the reaction rr + K + — > n + K + . This is studied 
explicitly at T = in Sec. IV A and the scattering lengths for this process as well 
as the remaining 7r.fr processes, which can be obtained via crossing symmetry, are 
then evaluated. In Sec. |1V H , the generalization to finite temperature is discussed. 
We summarize and conclude in Sec. [V|. 



II. SU(3) NAMBU JONA-LASINIO MODEL 

The SU(3) NJL model is well documented elsewhere ||||. In this section, we 
give only a brief description of the model in order to introduce our notation. Our 
starting point is the minimal three flavor chiral Lagrangian density 



^chiral = & W + G £ [0M» 2 + (^7 5 A» S 
a=0 

- K [det$(l + 76 )V) + det$(l - 7b )V>)] (2.1) 



that contains a four point coupling with strength G, and six point coupling, with 
coupling strength K. Color, flavor and spinor indices of the quark fields ip are 
suppressed. The term moderated by G contains the Gell-Mann matrices A a , a = 1..8 



plus the matrix A = J 2/31, where I is the unit matrix, thereby ensuring the 
U(3) symmetry of this term. This symmetry is broken by the determinantal term, 
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motivated by instanton effects, down to the chiral group SUl(3) x SUr(3). In order 
to incorporate the known current quark masses, an additional term that explicitly 
breaks the chiral symmetry, 



{m° u uu + m° d dd + m°ss) 



(2.2) 



must be introduced, so that £njl = £ c hirai + £ ma ss- 

A straightforward evaluation of the self-energy in the mean field approximation 
yields constituent quark masses that are modified from their current values due to 



dynamical symmetry breaking [19 



mi = m \ + AiGNctxS* - 2KN*(trS j )(trS k 



(2.3) 



where the propagator S l (x, x') satisfies 



(i0 x -m i )S i (x,x')=5 (4 \x-x f ) 



(2.4) 



in the mean field approximation and i is a flavor index. 

The mesonic sector is viewed most simply from the effective NJL Lagrangian 
that is obtained on constructing an effective four point interaction from the six 
point term, on contraction of two fermionic fields. One has 

8 



+ 
+ 



i=0 
1 



ijfc>$Aty)(W) + ^i +) (^7 5 AV)(^7 5 AV) 



(2.5) 



with the effective couplings 



(±) 



G =F -N c K{itrS s + 2ttrS u ) 



r[ ±] = Ki ±] = Ki ±] = G± ^N c KitiS s 

K[ ±] = = = K\ ±] = G± -N c KitrS u 

= Gt -N c K(itrS s - 4itrS u ) 
6 

K£ ] = t\n c K(HtS s - itiS u ) , (2.6) 

in the case m° = m° d . Both scalar and pseudoscalar octet channels may be identified 
by the isospin operator correspondence 

A 3 — > 7T°, a n 
-±=(\ 6 ±t\ 7 )^K ,K ,a K ,a K 

-L(\ A±l \ 5 )^K^,a£ ) , (2.7) 

using an obvious notation for the scalar sector. Experimentally, one could tentatively 
assign the scalar particles (0^,0^) — > (ao(980), -K"q(1430)). The assignment of the 



(er, a/) to actual particles such as the controversial /o(400 — 1200) or a |22| or / (980) 
and /o(1300) is uncertain and we thus retain our model notation here. The pionic 
and kaonic scalar and pseudoscalar sectors can be treated independently now. The 
quark-antiquark scattering amplitude, evaluated in the random phase approximation 
is 

2K {±) 

M na = 7-^-- (2.8) 

' * 1 - 2if{ ±) lT]? 5 (fc 2 ) 1 ' 

for the pions, or 

2K {±) 

Mk,. k = tJ—- , (2.9) 

1 - 2Ki ±) n s « s (k 2 ) 

for the kaons, expressed in terms of the irreducible polarization 

- iU f ^ s 2 (k 2 ) = -2N C [ A^t T [iS^(p + q)T p > s iS f2 (q)T p > s ] . (2.10) 
J (27r) 4 



Here f\ and /2 refer to the quark flavors q = u,d and s. The indices P, S refer 
to pseudoscalar or scalar functions respectively, and T p = ry 5 , while V s = 1. The 
meson masses with the appropriate quantum numbers are obtained by solving 



1 - 2K[ ±) U qq s (k 2 ) = 



1 - 2K[ ±) U sq s (k 2 ) = 



p,s\ 



(2.11) 



for the pionic or kaonic masses respectively, while the associated couplings are given 

as 



Qhl^qq 



dk 2 



9M K sq 



dk 2 



(2.12) 



In this equation, M w refers to the mesons n or a n while Mjf refers to the mesons K 
or ax as required. The mesons r] and 77', or correspondingly the scalars a and cr' are 
coupled. Here the quark-antiquark scattering amplitude in either sector takes the 
matrix form 



M = 2K(1-2UK)~ 1 



(2.13) 



where 



with 



n = 



/ n n \ 

WOO ^08 

nso n 8 8 



n 00 

nos = n 80 



^(2U qq + U ss ) 
\f2 

l—(U qq -u ss ) 
-(U qq + 2U SS ) 



(2.14) 



(2.15) 



and 
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K 



( K K \ 

-"■00 ^08 
-^80 -^88 



(2.16) 



Ku = Kf, K 08 = K 80 = X -K± (2.17) 

in either the scalar or pseudoscalar sector. For the purposes of our calculation of ttK 
scattering lengths, we will not require information on the r\ and rf mesons. However, 
the exchange of a a or a' in the intermediate state within itK scattering may occur. 
We thus require the quark-antiquark scattering amplitude for the exchange of the a 
and a'. This can be written as 

D ay = M 00 A x A + M 08 A x A 8 + M 80 A 8 x A + M 88 A 8 x A 8 . (2.18) 

Since the neutral a and a' are both exchanged in any single process, it is not nec- 
essary to resolve this scattering amplitude into the individual components for these 
particles via diagonalization. For our purposes, the individual masses and couplings 
are not required, although this provides an alternative method. The reader is re- 
ferred to Refs. PPP]|, who perform such a diagonalization explicitly. 



III. SCATTERING LENGTHS FOR ttK —> ttK AT T = 

In general, if one regards the process a + b — > a' + b', the scattering amplitude 
for a fixed isospin I can be decomposed into partial wave amplitudes, i.e. 

oo 

T 7 (s,t,w) = 16vr^(2/ + l)*f(s)i}(cos0) (3.1) 

1=0 

where s, t and u are the usual Mandelstam variables s = (p a + Pb) 2 , t = (p a — p a >) 2 
and u = (p a — Pb) 2 -, an d 9 is the scattering angle. The phases Si(s) parametrize the 
scattering amplitude as 

tf(«) = #^ \e 2l& '^ - l] , (3.2) 



2q 2i 
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where 



1 



q = ^~y|V I s ~~ ( m « + m b) 2 ] I s ~ ( m a - m b) 2 } ( 3 - 3 ) 

is the center of mass momentum of the incoming particles. The differential cross 
section is 

<h ' ' Vis, t) | 2 , (3.4) 



dil 64tt 2 s 1 

while the total cross section can be expressed purely in terms of the phase shifts, 

Air 

<** = 4 £( 2 ' + 1) sin2 ti(8) . (3.5) 
1 i 

At low energies, the partial wave amplitudes can be expanded in the form 



4(8)] =^<t»[aj + $<? + ()(<?)] , (3.6) 



defining the scattering amplitudes af and the slope parameters bf. At the kinematic 
threshold q — > 0, the cross sections in each isospin channel are given purely in terms 
of the s— wave amplitudes, i.e. 

a 1 = 4tt (off . (3.7) 



Our interest is in the nK system, and we start with the process of maximal 
isospin, I = 3/2, 

7T+ + K + -> 7T+ + K + . (3.8) 

By exchanging particles a and a', one directly obtains the scattering process 

n- + K + ^ it' + K + , (3.9) 

in which u is the center of mass variable. This reaction contains not only an I = 3/2 
component, but the 1 = 1/2 component also. Crossing matrices relate the various 
amplitudes. Using an obvious notation, one may write 

11 



3/2 

T'(s,t,u)= £ C^T^s^u) 

I'=l/2 

= C s I u 1/ X /2 (s,t,u) + C s I u 3/2 T^(s,t,u) , (3.10) 



with the crossing matrix |2_ 



Cl > r ~ 3 



(3.11) 



V 2 l J 



In practice, this tells us that 



TV\s, t, u) = 3 -T^(u, t, s) - V 2 ( S , t, u) (3.12) 

so that the determination of scattering in the isospin channel / = 3/2 is sufficient 
to obtain all information, including the amplitude in the T = 1/2 channel. We 
therefore examine the process ir + K + — > ir + K + in what follows. 

A. General Classification of Diagrams 

Since the coupling strengths in the NJL model, GA 2 and KA 5 are large, a pertur- 
bative expansion in the coupling strengths is inadmissable. An alternative expan- 
sion in the inverse number of colors N c is generally made |T0|JTl| . In this expansion, 
fermion loops contribute a factor of A^ c , while the scattering amplitudes with mesonic 
intermediate states contribute a factor of 1/N C . Thus the diagrams that are leading 
in the 1/N C expansion, and which are of the same order in this expansion, are both 
the box diagrams of the type displayed in Fig. [I] and the mesonic exchange graphs 
of Fig. g. 

B. Contributions to ir + K + — » ir + K + . 



12 



1. Box diagrams 



For the process tc + K + — > it + K + , only one of the box diagrams of Fig. [I] that is 
commensurate with isospin and strangeness conservation survives. This is depicted, 
togther with the appropriate kinematic variables, in Fig. |3|. A direct translation of 
this diagram leads to the expression 

iT hox = -4,{ig^ qq ) 2 (igKsq) 2 N c (3.13) 

/d q 
(2^)4 tr i ir y5iS 9 (q)i-f 5 iS q (q - Pi)i^iS q (q - k 2 + ki)ij 5 iS s (q - k 2 )] ■ 

Here the factor 4 arises from the flavor algebra. It has also been assumed that the 
translation of the Feynman graph gives %T . 



2. Meson exchange graphs 

As in the previous subsection, it is necessary to analyse first which diagrams can 
contribute to the ir + K + elastic scattering amplitude. One notes immediately that 
the s channel graphs of Fig. ^| do not contribute, due to charge conservation - the 
intermediate state would require a doubly charged (scalar) meson, which does not 
exist within the model. In the t channel, on the other hand, there are six possibilities. 
These are displayed explicitly in Fig. |j. Note that in all these diagrams, a scalar 
meson is exchanged in the intermediate channel. This is a consequence of parity 
conservation and manifests itself directly in the evaluation of the triangle diagrams 
that form a part of each amplitude in Fig. ^, in that the trace of the product 
of 75 675 675 S vanishes, whereas that of 755755'IS' does not. We construct these 
diagrams by first considering the triangle graph that is given in Fig. |5|. From this 
Feynman diagram, the three meson vertex is given as 

- iT^ih, k 2 ) = —N c J -^tr[i l5 iS\q)i lf> iS 2 {q - kJiS^q - k 2 )] (3.14) 
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fti + Tn 2 )i(ii— % + Tn 2 ) 



(2tt) 4 (g3 _ m 2) [(g _ fcl )2 _ m 2] [(g _ ^2 _ m 2] > 

where the superscripts 1 and 2 denote the flavors of the internal quarks and no trace 
over flavor has been included at this stage. This must be combined with the effective 
propagators of the scalar mesonic sector. We may simply use the form 



Da^ifa - k 2 ) = MqoAq x A + M 08 A x A 8 + M 80 A 8 x A + M 88 A 8 x A 5 



(3.15) 



where the My are a function of momenta, My = My(fci — k 2 ). If, for example, a 
strange quark couples to the first vertex and a u quark to the second of the mesonic 
amplitude, as is the case in (c) and (d) of Fig. |], one can obtain the flavor factors 
by considering the flavor components of Eq. (|3.15|) separately. The 8-0 term of this 
expression, for example, is thus 









M 


Xa{s,s) ■ X (u,u) = (001)A 8 





■ (100)A„ 














2V2 



and the other terms follow similarly, leading to the expression 
D 9 ' 



2 y/2 2 

-Moo - ~yM 08 - -M 88 



(3.16) 



(3.17) 



This form of the intermediate mesonic scattering amplitude is required for the eval- 
uation of diagrams (c) and (d). With it, one easily constructs the t channel txK 
scattering amplitudes in (c) and (d) to be 

i T t!{(c)+(d)}^ fc 2,Pi,p 2 ) = -i2 ■ 4c4 ? #iUrf Oi, h)D s a q a'( k i ~ k 2 )Vf{-p 2 , - Pl ) , 

(3.18) 

since terms (c) and (d) are equal when m u = m^. In this expression, the factor 4 
arises from the external vertices, that each contribute a factor y/2. The factor 2 
comes about since the contributions from graphs (c) and (d) are equal. 
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In an analogous fashion, diagrams (a) and (b) can also be evaluated. The inter- 
mediate mesonic scattering amplitude in this case is found to be 

DZ'(h-k 2 ) = ^M 00 + ^M 08 + ^M 88 , (3.19) 

and in combination with Eq. (|3.14j) leads to the expression 
*'^ 3 {(a)+(5)}( A; i ) h,Pi,P2) = -*2 • Agl qq g 2 Kqs Y\\-k 2l -k^D^fa - k 2 )Yf{p u p 2 ) 

(3.20) 

for the t channel tt + K + scattering amplitudes depicted in (a) and (b). 

There is no contribution from diagrams (e) and (f): the only difference between 
these two graphs is seen to be due to the flavor couplings to the intermediate ex- 
changed meson a n . Since \s(u,u) = —X 3 (d,d), the sum of these two terms cancel, 
and hence 

<(e) + (/)} =0 • ( 3 - 21 ) 

Thus, in the t channel, the tt + K + amplitude leads to 

rp3/2 rp3/2 rp3/2 /r, r\r\\ 

2 t - I tUa)+(b)\ + I tUc)+(d)} > 



With T tUa)+lb)\ aIld T t{(c)+(d)\ aS S iven 111 E( 1 S - O) aIld (EID 



t{{a)+(b)} a ^ J t{(c)+(d)} 



We now come to the u channel graphs of Fig. |2|. Such processes require the 
exchange of an uncharged strange meson. In the NJL model, this corresponds to 
the exchange of a <jk and can be realized within the model as depicted in Fig. 
An elementary vertex needed for this process is shown in Fig. [7]. Analytically, it can 
be constructed as 

-tTf{h 1 ,p 2 ) = -N e J ^triiKiS^q^iS'iq-kJiS^q-p,)] (3.23) 

_/V f d4q t t^M. + m i)^5%- h + m 2 )i{4- j> 2 + mi) 
c J (2vr) 4 1 (g 2 - m\) [(q - - m 2 f] [{q - p 2 f + m\] 
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while the effective interaction mediated by the ok meson follows from Eq. (2.S) as 



D (JK {k l -p 2 ) = r , (3.24) 

1 - 2K 44 Yl b sq (k 1 -p 2 ) 

so that one arrives at the form 

iT^ih, k 2 , Pl ,p 2 ) = -Sigl q y Ksq Tr(k 1 ,p 2 )D aK (k 1 -p 2 )rr(k 2 , Pl ) (3.25) 

for the u channel amplitude. The factor 8(= y/2 6 ) in this expression comes from 
all flavor factors that contribute at each vertex. 

The complete 7 = 3/2 tt + K + — ► tt + K + scattering amplitude is thus made up of 
three components, 

T 3/2 (k l ,k 2 ,p 1 ,p 2 ) = T^(ki,k 2 ,p l ,p 2 ) +T^ /2 (k 1 ,k 2 ,p 1 ,p 2 ) + T 3/2 (k l ,k 2 ,p l ,p 2 ) , 

(3.26) 



with T box , T t 3/2 and T^ 2 given by Eqs. (gTOf ), ( ggg ) and ( |3725| ) respectively. 



C. 7TK scattering lengths 

The expressions given in the previous section enable one in principle to calculate 
the 7rK cross section for an arbitrary choice of kinematic variables. In practice, this 
turns out to be extremely difficult, as the box diagram cannot easily be evaluated 
exactly for arbitrary kinematics, requiring in general addition approximations. In 
this paper, we restrict ourselves to a calculation of the scattering lengths, which 
simplifies the calculation somewhat. The kinematic threshold is given by 

s = (m 7r + m^) 2 (3.27) 

with 

t = 0, u={m n -m K ) 2 , q = . (3.28) 
16 



These conditions can be fulfilled by choosing 



k\ — &2 = k 



m K 




(3.29) 



and 



Pi = P2 = P 







(3.30) 



/ 



Using these kinematics, the box term, as well as the triangle and intermediate meson 
exchange graphs making up the t and u channel contributions, can be expressed in 
terms of the "elementary" integrals 

d 4 q 1 



F' 
M 12 (p) 
N 12 (p) 
P 12 (p,k) 
Q 12 (p,k) 



(2tt) 4 q 2 - ml 
d 4 q 



1 



(2tt) 4 (q 2 - ml)[(q - p) 2 - m|] 

d A q 1 
(27r) 4 (q 2 — mi) 2 [(q — p) 2 — m?,] 

d A q 1 
(27r) 4 (q 2 — mf) [(q — p) 2 — m 2 ] \(q — k) 2 — m 2 ] 

d 4 q 1 
(27r) 4 (q 2 — m 2 ) 2 [(q — p) 2 — m 2 ] [(q — k) 2 — m\ 



(3.31) 
(3.32) 
(3.33) 
(3.34) 
(3.35) 



that form the building blocks for these functions. We list the results for the compo- 
nents of the T 3 / 2 amplitude for the threshold kinematics. 



1. Box diagram 



The kinematic choice of Eqs. ( |3.29| ) and (|3.30|) lead to the forms 



iT^(k,k,p,p) = -AN c gl qq g 2 Kqs 



(3.36) 



x 



/ 



d 4 q / tr [7 5 (jf- j£ + m s )-/ 5 (d + m q )^{$[- p 1 + m q )^(4 + m q )\ 



[(q — k) 2 — m 2 ] [q 2 



m 



212 [(q — p) 2 — m 2 } 
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After some calculation along the lines that will be indicated for the three meson 
vertex function to follow, one finds 



~ 8gl qq 9 2 Kqs N c {M»(0) + M" s {k -p)- p 2 N™(p) 

- ~k 2 - (m q - m s ) 2 ] N qs (k) - 2k ■ pP qs (p, k) 

+ p 2 [k 2 - (m, - m s ) 2 ] Q qs (p, k)} . (3.37) 



2. Meson exchange graphs 

The three meson vertex functions are required for the t and u channel graphs. 
From Eq. (|3.14|), one has 



1V ' ; J (2vr) 4 [q 2 - m 2 ] 2 [(q + k) 2 - m 2 } V ; 

which reduces to 

- iT l2 (k k) - -UN f A g 2 (^i ~ 2m 2 ) - 2m 2 q ■ k + m x m 2 

after performing the trace. Using the relation 

2k ■ q = [{k + q) 2 - m\] - [q 2 - m 2 .] - k 2 -m\ + m\ (3.40) 

leads one to the expression 

- iTf{k, k) = -AiN c [ - m 2 M 22 (0) + (mi - m 2 )M 21 (k) 

+m 2 [k 2 - (m 2 - mi ) 2 }N 21 (k)} (3.41) 

in terms of the elementary integrals. In a similar fashion, the second vertex function 
that is required for the u channel T matrix amplitude of Eq. fl3.23j ) and which was 
shown in Fig. |7| can be decomposed as 
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iT 12 (k,p) = 2N c i[(m 2 - mi)M 12 (A;) + (m 2 + mi )M 12 (k - p) 
-[m 2 p 2 + m l k 2 -m l (p-k) 2 ]P 12 (p,k)} . 



(3.42) 



In addition to this, the scalar polarization function is given as 
- iU 12 (k) = -2N J + ft + m2 ^ 

and is expressed as 



(2tt) 4 [q 2 -m 2 ][(q-ky -ml] ' 



(3.43) 



ill 12 (k) = 4iV c {F 1 + F 2 + [(mi + m 2 ) 2 - k 2 ]M 12 (k)} . 



(3.44) 



Note that the pseudoscalar polarization differs from this function only in the re- 
placement (mi + m 2 ) 2 — * (mi — m 2 ) 2 . 

The complete expression for the I = 3/2 scattering amplitude at the kinematic 
threshold follows now on constructing 



ilfiL , fM (k, k,p,p) = -d>igl qq g 2 Tf{k, k) I ^M 00 + ^M 08 + ^M 88 ] rf (p,p) 



2^2 



t{(a)+(b)y 



(3.45) 



and 



V2 



rf(p,p) 

(3.46) 



2fG 



for the t channel, and 

iT* /2 (k, k,p,p) = -8igl qq g 2 Kqs [If (k,p)} 2 1 _ 2 ^nf 8 (fc-p) 
for the u channel, with the couplings 



(3.47) 



gj c i a = -4iNAM»(k) + 



(m q — m s ) 2 — k 2 



2k 2 

+ (m 2 s -m 2 + k 2 )N sq {k) 



M qs {k) - M ss (0) 



(3.48) 



k 2 =m 2 K 
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and 



g~ 2 qq = -2iN c 



M w (0) + M qq (k) - k 2 N qq (k) 
which can be derived from the definitions in Eq. ( p. 12 ). Then one has 



T 3/2 (k, k,p,p) = T^ 2 (k, k,p,p) + Tf/? )+m (k, k,p,p) 



+ T t{(c)+(d)} fc > P> P) + T u /2 ( k , k , V, V) ■ 



The scattering lengths are evaluated as 



T 7 (s = (m n + m K ) 2 , t = 0, u = (m n — m K ) 



(3.49) 



(3.50) 



(3.51) 



8^(771^ + 777^) 

so that a 3//2 follows directly from Eq. ( |3.5U| ), while a 1 ' 2 is obtained subsequently 
through the relation Eq. (|3.12|). 



IV. NUMERICAL EVALUATION OF THE nK SCATTERING LENGTHS 



A. T=0 Case 

In order to evaluate the scattering lengths, the functions F, M, N, P and Q of 
Eqs. (|3.31|) - (|3.35| ) must be evaluated. Since the NJL model is non-renormalizable, 
a cutoff parameter must be introduced to regulate these functions. This can be 
done in a covariant fashion using the Pauli-Villars scheme, or non-covariantly by 
introducing a cutoff on the three momentum |p] < A3. For consistency, a cutoff 
is included in all integrals, including N, P and Q, that are convergent. Technical 
details for the Pauli-Villars integrals are given in Appendix [A[ The 0(3) calculation 
is performed directly numerically. Our first evaluation is performed at T = 0. 

In Table H, we list two sets of parameters for the Pauli-Villars calculation that 



were taken from Ref. [|T9 | , as well as two sets of parameters for an 0(3) cutoff, taken 
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from Refs. @,^(J. These calculations are all non-chiral, all assuming non- vanishing 
values for the current quark masses. Using these parameter sets, we obtain values for 
the pseudoscalar meson sector at T = that are listed in Table [II]. The experimental 
values are also given. The scalar meson masses for T = for the corresponding 



parameter sets are given in Table |T|. Using these values, the s-wave scattering 
lengths a]J 2 and (Xq £1X6 calculated at T = 0. The results for a^ 2 , obtained from 
the box, t and u channels for each parameter set, are given in Table [TVJ with those 
for a]j 2 given in Table 0. The final results for the scattering lengths arise from a 
delicate cancellation of the box, t and u channel values, which is evident from these 
tables^ In the model, we find that the calculated values of a^ 2 range from —0.003 

1/2 

to —0.04, while that of a range from 0.13 to 0.18, in units of m n . Our results 

1/2 

for a fall well within the experimentally known range of 0.13 to 0.24, while those 
for a^ 2 are lower, the experimental values being —0.13 to —0.05. Standard chiral 
perturbation theory gives a^ 2 = 0.17 and a^ 2 = —0.05 [|TJ for these values. Thus 
our T = values, while dependent on the regularization scheme, can be considered 
reasonable. 

B. ttK scattering lengths for T / 

At finite values of the temperature, the irK scattering graphs can be analysed 
as was done at T = 0, but within the imaginary time or Matsubara formalism. This 
corresponds to making the formal substitution 



3 To obtain a deeper understanding of this fact, it is necessary to perform a chiral ex- 
pansion of the box, s, t and u channel graphs in the two variables (m n /A and ra^/A) 



analytically along the lines suggested by Refs |17| and [|18[. This is a difficult problem in 
its own right and will not be addressed here. 

21 



dtq "zf.^ («: 



(27,)' /} V Ja (2t) 3 
and replacing the fermionic (bosonic) go by the discrete frequencies zu; n , co> n = (2n + 
1)tt//3 or u; n = 2nn/i3, n = 0, ±1, ±2, ±3 • ■ ■ respectively. The final expressions 
are given in Appendix B. As is indicated in Eq. (|4.1| ), an 0(3) cutoff is naturally 
implemented. 

The pseudoscalar meson pole masses vary with temperature, and this variation 
is shown in Fig. |[ together with the dissociation threshold for the kaon, m q + m s , 
which leads to a kaon dissociation temperature Tmk- The function tuk + Tn n is 
also shown. The temperature at which this function crosses m q + m s defines the 
point T n x at which a meson pair can dissociate into a quark pair. Numerically, for 
the set of parameters denoted as HK in Table 0, we have T n x = 174.5 MeV, while 
Tmk = 203.0 MeV. We also show the temperature dependence of the scalar mesonic 
sector in Fig. |9] for this parameter set. Note that the behavior of the a n meson 



with temperature p0| has been confirmed qualitatively in lattice gauge calculations 
|13| . Using the temperature dependent functions for the mesons, the temperature 
dependence of the box, t and u channel amplitudes is evaluated and used to construct 
the scattering lengths. In Fig. [K], we show the individual contributions from these 
channels, together with the evaluated scattering length a . One sees that all 
graphs diverge close to the Mott point. The box and u channel graphs do not 
change substantially up until the Mott point, while the t channel meson exchange 
diagram displays a somewhat faster variation with temperature. Consequently, the 
resulting change in a^ 2 is noticeable - the change, within this model calculation, 
is quantitatively given as a^^T = 150 MeV)/a^ 2 (T = 0) = 2.36, in contrast to 



the completely flat temperature dependence found for nn scattering lengths pi 
The temperature dependence of the a^ 2 is less dramatic: Fig. shows that the 
sum of the functional components lead to a flat behavior. Quantitatively, one has 
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a /2 (T = 150 MeV)/a /2 (T = 0) = 1.18. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have evaluated the temperature dependence of the 7iK scatter- 
ing lengths a^ 2 and clq 2 within the framework of the SU(3) NJL model. Ordering 
terms according to an expansion in 1/N C is used to determine the leading set of 
diagrams that contribute to the scattering amplitude. One finds a box-like graph, 
as well as t and u channel graphs. To this order in N c , no unitarity corrections are 
present. Because of this, only the scattering lengths are investigated, although the 
formulae are given for arbitrary kinematics, so that the scattering amplitude can 
be evaluated for arbitrary kinematics. The restriction here to threshold kinematics 
also simplifies the computation considerably. 

At T = 0, we have calculated the / = scattering lengths a^ 2 and a?J 2 for 
different regularization schemes, that of Pauli and Villars, and also using an 0(3) 
cutoff. We have simply used parameters sets of other authors and have not fitted 

1/2 

any additional quantities. We find that a falls well within the experimentally 
known data bracket, while a^ 2 is slightly lower than the bracketing values given 
experimentally. Unfortunately the experimental values are not known precisely at 
this stage. 

In investigating the finite temperature behavior of the scattering amplitudes, we 
observe that the box and u channel graphs, the latter of which occurs by resonance 
exchange of a a K scalar meson, do not vary strongly with the temperature. However, 
variation of the t channel amplitude that is controlled by the exchange of o and a 1 
mesons displays a slightly more pronounced temperature dependence, leading to a 
temperature change of a^J 2 of a factor of about 2 over 150 MeV. This is in contrast 
to the case for rnr scattering, in which the hit scattering lengths are found to be 
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essentially independent of temperature over a range of about 150 MeV. The study 
of kaon physics appears here to give a stronger functional dependence than in the 
pionic case. Since this is driven essentially by the temperature dependence of the 
intermediate scalar mesons, it seems that processes that involve the exchange of 
mesons in the scalar sector may ultimately provide us with a clean signal for a 
chiral phase transition. 
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APPENDIX A: PAULI-VILLARS REGULARIZATION 



Standard Feynman parametrization ||24j| , 

17— = SO f 1 ffr dx . #A *(!-*) 



(Al) 



1 ) ixUx t A t Y ' 

with ( = Yli=id an d x = J27=i x i i s used for the integrals Eqs. ( |3.31| )- fl3.35| ). All 
calculations are performed at the kinematic threshold. One has 

d A q 1 



(2vr) 4 q 2 - m 2 ' 



(A2) 



which, after Wick rotation, gives 



(4tt) 



2 2 i I 2 2 

q — m m\q + m 



CO 





(A3) 



The Pauli-Villars regularization scheme is then implemented: 
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■— - ]T Cj [(m 2 + a 3 A 2 ) In \q 2 + m 2 + a,A 2 







im 



(Any 



In 



m 

J 2 



-2 



+ 1 



A 2 
2— In 



( m 2 \ 1 

1+ U +1 



(A4) 



where m 3 - = m 2 + ajA 2 and the standard set of parameters Co = 1, C\ — 1, Ci — — 2, 
«o = 0, a.\ = 2 and a2 — 1 are used. 

The function M 12 (p) from Eq. ( |3.32| ) can be integrated after Wick rotation. One 

has 



M l \p) = i 



d 4 q E 



dx 



(2tt) 4 Jo (q 2 + R 2 f 



(4vr) 2 Jo 



dx 



R 2 



q 2 + R 2 



ln(g 2 + R 2 



(A5) 



J o 



with R 2 = m\x + m 2 (l — x) — p 2 (l — x)x. The Pauli-Villars regularized version reads 



M 12 (p) 



W?]k Cj h dx 



R2 > + \n(q 2 + R 2 ) 



_q 2 + R) 



(A6) 



J o 



with the replacements m\ — > m\ 3 = m\ + <x,A 2 and m 2 —>■ m\i = m\ + <x,A 2 in R 2 



defining R 2 . One obtains 



M 12 (p) 



i 2 r 1 
— Cj / dx In 
lvr 2 £^ 3 Jo 



(m\ — m 2 )x + m 2 j — p 2 (l — x)x 



The ^-integration can be performed explicitly leading to 



(A7) 



M 12 (p) 



(4tt) 2 



i=o 



2p 2 



In 



m 



1.7 



2j/ P 



[At 



with Aj = (m\ — — p 2 ) 2 — Am\ 3 p 2 and 
A,- = 



-A, 



tan" 



-1 / P 2 + m l — m 2 



tan' 



-A, 



(A9) 



if Aj < and 



A i = \/ A i 



1 - my - m 2 ^ + JAj 
- In 5L__ 

2 ^2 _ m 2^ _ m 2^ _ / ^. 



inQ (p 2 — (mi — m 2 )' 



(A10) 



25 



if Aj > 0. Following the same prescription for N 12 (p) leads to 



dx 



x 



(4tt) 2 Jo R 2 



(An) 2 2p 2 



ln(^]+2 



mj — — p 



A 



rrio 



(All) 



Both P 12 (p, k) and Q l2 (p, k) of Eqs. fl3.34| ) and ( |3.35| ) are functions of two vari- 
ables and are somewhat more complicated. Feynman parametrization followed by a 
Wick rotation and integration over the momentum variable leads to 



(4vr) 2 Jo 



i ,i 

dx / ay 
Jo 



y_ 

R 2 



(A12) 



where 



R = m 1 (l — y + xy) + m 2 (l — x)y — [k (1 — x) +p x]y(l — y) — x{\ — x)y u 



(A13) 



with u = (p — k) 2 . The integral in x can again be performed. One finds 



P u (p,k) 



2i 



dy 



tan- 1 f 2u y+p_-*y I _ tan -i ip-«y 



(An) 2 Jo dA 
where the discriminant in this case is 

- y 2 & = (pV - ny 2 ) 2 - 4(m 2 + ujy + k 2 y 2 )y 2 u 



■ (A14) 



-y 



(4m 2 u - p 2 ) + (2pn + Auu)y - (k 2 - Ak 2 u)y 2 



(A15) 



with p = m\ — mj + k 2 — p 2 , k = k 2 + u — p 2 and to = m 2 — m\ — k 2 . 
A similar analysis leads to 



Q U (p,k) 



dxdy^ 



(A16) 



(An) 2 Jo Jo (fry 

in which the x integration can once more be carried out with the aid of the integral 

_! / 2ax + 6 N 1 ' 



dx 1 
o [ax 2 + bx + c] 2 A 



2ax + b 4a 

— 7. — ; 1 — ; tan 

ax 1 + bx + c v — A 
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(A17) 



with A = b 2 — Aac. Note that the Pauli-Villars regularization procedure has not been 
carried out explicitly on N 12 (p), P l2 (k,p) and Q 12 (k,p), which are convergent. For 
consistency within the model, however, this should be done. The implementation 
is simple. The standard replacement of each integral by a sum of three terms with 
Cj and rriji is constructed and evaluated numerically. The forms for P and Q, 
Eqs.(A14) and (|A17|) are valid in the case that the integrands are non-singular. 



APPENDIX B: FINITE TEMPERATURE INTEGRALS 



Using the abbreviation E 2 = q 2 + m 2 for the energy of particle with flavor j, 
one finds the following forms for the integrals Eqs. fl3.31p - ( |3.35| ) at the kinematic 
threshold: 



F j 



i f A , q 2 



4tt 2 Jo * Ej 



dq— tanh 



while 



M 12 (p) 



167T 2 Pq 



Qo - ¥ 



tanh | 
+ tanh 



/3Ei\ pl + m\ — m\ 

~T) E x 

/ (3E 2 \pl + m 2 2 -ml 



—i- 



nq 



tanh 



+ tanh 



v 2 ; E 2 

' PE W \ gl + m 2 - m 2 
, 2 j E 10 
'f3E 20 \ ql + m 2 2 -m\ 



\ 2 J E 2 o 
where we have abbreviated 

[pi - (mi + m 2 ) 2 ][p 2 - (mi - m 2 ) 2 ] 



4 



Apl 



E 



jo 



ll + m 2 



(Bl) 



fB2l 



(B3) 
(B4) 



For N (p) we find 
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i 1 f A , q 1 
dq- 



32n 2 p%Jo ^(q 2 -q 2 f 



x 



tanh I 
— tanh 



2 / 4^ 
'f3E 2 \ (p 2 -m 2 + m 2 2 ) 2 + 4p 2 E 2 



2Eo 



13 Po {% ~ (?) (Po + m l ~ m l 



cosh 2 ( 



f3E 1 



2E{ 



(B5) 



with £ = p\ + to 2 — to 2 ,. For the function P 12 {k,p) one has 



P 12 (p,A;) = 



<ig- 



47r 2 7o 5 p 5fc5, 



tanhf^ ^ + ^ + tanhf^^ 
2 E 1 \ 2 E 2 



(B6) 



with 



5 P = 
4 = 
S u = 
A = 
B = 
C = 



(A; 2 + to 2 - to 2 ) 2 - 4£ 2 £; 2 = (A; 2 - to 2 + to 2 ) 2 - 4£ 2 fc 2 1 



■u + m\ — to 2 ) — AE 2 u 
pi (kl + to 2 - m 2 ) + 2E\{p\ + k 2 -u) 
pi (u + m\- to 2 ) + 2E\{p\ - kl + 



M — TO 2 + TO 2 ) — AE\l 



(k" 2 , — to 2 + to 2 ) (w — to 2 + m|) — 2E 2 (pl — Ajq — -u) 
For the function Q 12 (k,p), one has 



(B7) 
(B8) 
(B9) 
(BIO) 
(Bll) 
(B12) 



Q 12 (p,k) 



An 2 



— tanh 



2 J E 2 5 2 5 u 4 cosh 2 (fi)^A. 



(B13) 



where 



I? = p 2 [77 (2El + p 2 ) + 2£ 2 (4p 2 + 3A; 2 - 3m 

x '{pi + AE 2 ) ( V 2 + 4A; 2 £ 2 ) + 8 V E 2 (p 2 Q + fcg - «) 
- 4£ 2 [2^77 + (p 2 + A; 2 - «) (P 2 + pi) + 2p 2 £ 2 
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and 



£ 



2pl (r, 2 



Ak 2 E 2 ) + V (pl + k 2 - u) (pg + AE\ 
(pi + AEl) (u + m\- mj) + AE\ (p 2 - k\ + u) 
(k% -m\ + ml) 2 (u-m\ + m 2 2 ) + Ak 2 E\ (u-m\ + 
AE 2 2 (k 2 -ml + m 2 2 ) (p 2 - k 2 - u) 
2E\ (pi + k 2 -u)+p 2 (k 2 + m\- ml) 



m. 



(B14) 
(B15) 

(B16) 
(B17) 



; 2 i 2 2 

rj = k + m 1 — m 2 



(B18) 



These expressions simplify somewhat in the case mi = m 2 . 
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TABLES 

TABLE I. Parameter sets. The parameter sets denoted I and II correspond to those 
from Ref. jn| and which are similarly denoted therein. PR and HK designate the 0(3) 
parameter sets used in Ref. and Ref. || respectively. 





Renormalization 


m 0q (MeV) 


m 0s (MeV) 


A (MeV) 


OA 2 


KA 5 


I 


Pauli-Villars 


6.2 


175 


795 


2.350 


27.83 


II 


Pauli-Villars 


7.8 


175 


700 


2.739 


43.28 


PR 


0(3) 


5.5 


140.7 


602.3 


1.835 


12.36 


HK 


0(3) 


5.5 


135.7 


631.4 


1.833 


9.288 



TABLE II. Quark and pseudoscalar meson masses at T = 0. The parameter set labels 
I, II, PR an d HK are as described in Table |. Experimental values are given in the final 
row. 





m q (MeV) 


m s (MeV) 


m n (MeV) 


m K (MeV) 


m v (MeV) 


(MeV) 


I 


253.0 


489.7 


140.8 


521.9 


467.6 


845.7 


II 


367.8 


557.4 


140.0 


485.9 


497.3 


1013.4 


PR 


367.7 


549.5 


135.0 


497.6 


514.8 


957.7 


HK 


334.7 


527.4 


138.0 


495.7 


486.7 


857.4 


Exp 






138 


495 


548 


958 
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TABLE III. Scalar meson masses at T = 0. The parameter set labels I, II, PR and 
HK are as described in Table |. A tentative assignment of these masses to experimentally 
measured particles is given in the final row. 





m „ (MeV) 


m aK (MeV) 


m a (MeV) 


77v (MeV) 


I 


688.8 


904.6 


510.7 


1112.2 


II 


890.8 


1073.4 


735.4 


1226.3 


PR 


880.2 


1050.5 


728.9 


1198.3 


HK 


792.4 


980.0 


667.8 


1149.9 


Exp 


a (980) 


i^(1430) 


/ (980) 


/o(1300) 



TABLE IV. Contributions to the scattering length a in units of m n at T = 0. The 
parameter set labels I, II, PR and HK are as described in Table |. 



Parameter Set 


box 


t channel 


u channel 


sum 


I 


-0.3689 


0.2224 


0.1091 


-0.0375 


II 


-0.6454 


0.4577 


0.1587 


-0.0289 


PR 


-0.7047 


0.5332 


0.1741 


-0.0025 


HK 


-0.6000 


0.4000 


0.1771 


-0.0230 



TABLE V. Contributions to the scattering length a in units of at T = 0. The 
parameter set labels I, II, PR and HK are as described in Table |. 



Parameter Set 


box 


t channel 


u channel 


sum 


I 


-0.2255 


0.2224 


0.1296 


0.1265 


II 


-0.4706 


0.4577 


0.1447 


0.1318 


PR 


-0.5108 


0.5332 


0.1572 


0.1795 


HK 


-0.4195 


0.4000 


0.1710 


0.1515 
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FIGURES 

FIG. 1. Six possible box diagrams that contribute to meson-meson scattering. Double 
lines denote mesons, while the internal single lines depict the quark structure. 

FIG. 2. Three possible meson exchange graphs, corresponding to s, t and u channels 
respectively. The double lines represent mesons and the single lines represent quarks. 

FIG. 3. The only box diagram contributing to the elastic scattering amplitude 
n + K + — ► ir + K + . Here we denote the mesonic states with double lines and quarks by 
single lines. The s quark is denoted by a heavy single quark. 

FIG. 4. Possible t channel meson exchange diagrams that contribute to the scattering 
tt + K + — > tt + K + . Mesonic states are denoted by double lines, quarks by single lines. The 
s quark is denoted by a heavy single line. 

FIG. 5. Generic three meson vertex function T\ 2 that occurs in the t channel. The 
external double lines denote meson states. The internal single lines denote quarks with 
masses mi or mi as indicated. 

FIG. 6. u channel graph for the scattering amplitude tt + K + — > ir + k + . External double 
lines represent mesons, while single lines represent quarks. The s quark is represented by 
a heavy single line. 

FIG. 7. Generic three meson vertex function T^ 2 that occurs in the u channel. External 
double lines represent mesons, while single lines denote quarks with masses mi or ni2 as 
indicated. 
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FIG. 8. Pseudoscalar meson masses versus temperature for the parameter set labeled 
HK of Table |. Shown are m n (solid line), m,K (dashed line), and m v (lower dot-dashed 
line). The combinations m n + mx (dotted line) and m q + m s (upper dot-dashed line) are 
also shown. T n K and Tmk are indicated with arrows. 

FIG. 9. Scalar meson masses versus temperature for the parameter set labeled HK. 
Shown are m a (dotted line), m aiT (solid line), m aK (dashed line) and m a i (dot-dashed 
line). 

3 /2 

FIG. 10. Temperature dependence of well as the decomposition of the t, u and 

box graph contributions to this quantity for the parameter set of HK, in units of m n . The 
full curve for a^ 2 is indicated by the solid line, the box diagram contribution by dashed 
lines, the t and u contributions by a dotted and dot-dashed curve respectively. 

1/2 

FIG. 11. Temperature dependence of a as well as the t, u and box graph contribu- 
tions to this quantity for the parameter set of HK and in units of m n . The full curve for 

1/2 

a ' is indicated by the solid line, the box diagram contribution by dashed lines, the t and 
u contributions by a dotted and dot-dashed curve respectively. 
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